Underlying Event Studies at RHIC 



Underlying Event Studies at RHIC 

Helen Gaines for the STAR Collaboration 

Department of Physics, Yale University, New Haven, CT 06520, USA 



By studying p-p collisions we hope to improve our understanding of the fundamental constituents of matter 
and how they form into colorless objects. Measurements of the inclusive jet cross-sections and fragmentation 
properties have confirmed that QCD based calculations give a good description of the hard scattering processes. 
However, as our analysis of jets has improved it has become clear that there is significant contribution to 
these measurements from processes other than those directly related to the initial hard scattering - the so-called 
underlying event. Several processes contribute to the underlying event, namely the beam-beam remnants, initial 
and final state radiation and multiple parton interactions. The structure of the jet and the underlying event are 
strikingly different in both their particle compositions and momentum distributions. Only by understanding 
both components can one fully describe a p-p collision. I will discuss preliminary results from studies of the 
underlying event in p-p collisions at = 200 GeV at RHIC, and compare to PYTHIA predictions, as well as 
earlier results from the Tevatron at 1.8-1.96 TeV. 



1. Introduction 

The study of the properties of jets and the underly- 
ing event in p-p coUisions is important for improving 
our understanding of QCD and the hadronization pro- 
cess, as well as providing a vital baseline for compar- 
ison to measurements being performed at RHIC by 
the STAR and PHENIX collaborations in heavy-ion 
collisions [H El El S] . 

The results presented here are a preliminary study 
of p-p coUisions at = 200 GeV by the STAR col- 
laboration. It has previously been shown ^ that the 
inclusive jet cross-section in p-p collisions at ^/s = 
200 GeV is well described over seven orders of mag- 
nitude by NLO pQCD calculations. Further, these 
theoretical calculations can also successfully describe 
the TT and (anti)proton production in minimum bias 
collisions at the same collision energy [71 [8j. These 
results indicate that we have a well calibrated probe. 
However, the fragmentation functions, jet shapes, par- 
ticle composition of the hadronized fragments and the 
underlying event all remain unexplored at RHIC en- 
ergies. Such studies are now underway and a few se- 
lected preliminary results are presented below. 

Calculations suggest that the identified particle 
fragmentation functions are modified in different ways 
by the presence of a Quark Gluon Plasma [9 , this 
makes it especially important to have a well measured 
baseline with which to compare to the measurements 
in Au-Au collisions. 

These data are also compared to predictions of 
PYTHIA 6.410 [10], tuned to the CDF 1.8-1.96 TeV 
data (Tune A). This tuning has been shown to give 
good agreement with most of the measurements made 
at the Tevatron [TT]. An important test of whether 
such a tuning is truly representative of the underly- 
ing physics in elementary collisions, or merely a dili- 
gent selection of numerous parameters, is to compare 
PYTHIA predictions at a different collision energy 
without changing the parameters except via the es- 
timated energy scaling. A cartoon of many of the 



major components included in PYTHIA 's modeling 
of an elementary collision is shown in Fig. [Tj One 
such important energy scaling parameter in PYTHIA 
is the hard scattering cut-off of multiple parton inter- 
actions (MPI) in the collision. This cut-off is placed at 
PTo(Ec™) = PTo(Ecm/Eo)% where Eq is 1.8-1.96 TeV, 
the tuning collision energy. The initial estimate of e= 
0.16 has been widely used in LHC predictions, how- 
ever improved data at ^/s = 1.8-1.96 TeV and newer 
data at y/s — 630 GeV now suggest that a value of e— 
0.25 is more appropriate |llj . The inclusion of mul- 
tiple parton interactions increases the activity, and 
hence particle production, in the underlying event. In- 
creasing e creates a larger energy dependence of the 
multiple parton interactions. A change in PYTHIA 
of e from 0.16 to 0.25 results in a 35% increase in 
the charged particle density in the underlying event 
at RHIC energies, while decreasing the particle den- 
sity at the LHC by 26%. It is therefore important to 
determine the correct scaling parameter in order to 
have accurate predictions for the LHC when it begins 
collisions later this year. 



2. The Analysis 

We utilize data from the mid-rapidity Time Pro- 
jection Chamber (TPC) to reconstruct charged par- 
ticles, and the mid-rapidity Barrel Electromagnetic 
Calorimeter (BEMC) to measure the neutral parti- 
cle production. A "jet patch" trigger requiring Et >8 
GeV in a A77 X A0 = 1x1 patch of the BEMC was 
used to collect the data. During the RHIC Run-6 in 
2006 the use of this trigger allowed STAR to sample 
a luminosity of ^8.7 pb~^. By requiring such a large 
signal in the BEMC a neutral energy fragmentation 
bias is created for the triggering jet. Charged particle 
fragmentation functions are therefore presented only 
for the di-jet partner not associated with the triggered 
jet-patch. 
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Proton-Proton Hard Scattering Event 




Figure 1: Cartoon illustrating the major components used 
in QCD Monte-Carlo models to simulate a 2-to-2 parton 
hard scattering in a proton-proton collision. The event 
consists of two outgoing jets from the hard scattered par- 
tons and the underlying event. The jets contain particles 
from the fragmentation of the initial partons plus initial 
and final state radiation. The protons break up into beam- 
beam remnants. PYTHIA also includes the possibility 
of the second semi-hard 2-to-2 parton-parton scattering 
(multiple parton interactions). The underlying event con- 
sists of the beam-beam remnants, hadrons resulting from 
the multiple parton interaction and the initial and final 
state radiation hadrons not clustered in the hard scattered 
jets. Figure adapted from [5]. 



The kx and anti-ky and SISCone jet algorithms 
from the Fast Jet package fT^ were used to reconstruct 
jets, all three algorithms combine particles in rj — (j) 
space and are infrared and coUinear safe. SISCone, a 
seedless cone algorithm, merges particles within a cone 
radius R=y/A(j)'^ + Arj'^. A process of splitting and 
merging of cones is then performed to improve on the 
effectiveness of the cone algorithm in collecting energy 
radiated off of the scattered parton via higher-order 
QCD processes. The k^ and anti-k-r jet finders are se- 
quential recombination algorithms. Such algorithms 
successively cluster in a pair-wise fashion, all parti- 
cles and proto-jets (previously merged particle pairs) 
in the event via a pt and distance weighted cut selec- 
tion until no cluster pairs pass the cut. All remain- 
ing clusters are designated jets. For the k^ algorithm 
the selection starts from the lowest dy particle/cluster 
pair, where dij = m,in{k'^^, k'^j){A(j)^j + Ariij)^)/R^, 
R is the resolution parameter and k^ == Pt- If this 
minimum distance dij < /c^^, the pair are merged, 
and the iteration over all the new d^- continues un- 
til the test fails. The anti-k^ algorithm is similar 
except it merges from the highest dij and selects 
for = rmn{l/kl^,l/kj.^){A^l + Ar,fj)/R\ In 
this way the anti-k^ algorithm behaves much like an 
idealized cone algorithm. A cut of pt>0.2 GeV/c 
was applied to all charged particles considered in the 
event, and Et>0-2 GeV for each tower cluster re- 
constructed in the BEMC. To investigate how the 



radius/resolution parameters, R, used in this study 
aS'ect the reconstructed jet energy, the SISCone algo- 
rithm was first run with R=l. The energy contained 
within this jet cone as a function of R was then stud- 
ied. Fig. [2] it can be seen that higher energy jets are 
focussed within smaller jet radii. For jets with 20 
GeV/c < pi^^ < 30 GeV/c >75 (95)% of the jet's 
energy is contained within R=0.4 (0.7) , while >89 
(98)% of the jet's energy is contained within R=0.4 
(0.7) for 40 GeV/c < p!^* < 50 GeV/c. 




Radius of Jet Cone 

Figure 2: Color online: The fraction of the jet's en- 
ergy contained within jet cone radius R. Jets were initially 
found using the SISCone algorithm with R=l. 



The data are not yet corrected to the particle 
level, therefore the results are compared to PYTHIA 
6.410 calculations passed through STAR's detector 
simulation and reconstruction algorithms. The single 
charged particle reconstruction efficiency in the TPC 
is > 80% for Pt >1 GeV/c. Such detector inefficien- 
cies and the presence of particles which are undetected 
by the STAR detector, such as the neutron and K° , 
cause the reconstructed jet px to be lower on average 
than the true value. The errors shown in all plots are 
therefore statistical only at this point. 

The jet energy resolution was obtained via two tech- 
niques. The first uses the PYTHIA simulations to 
compare reconstructed jet energies at the particle and 
detector level. The second studies the energy balance 
of "back-to-back" di-jets in the real p-p data. Figure[3] 
shows that both methods resulted in comparable jet 
energy resolutions of ~20% for reconstructed jet pt 
> 10 GeV/c. 



3. Jets in p-p 

The uncorrected charged particle fragmentation 
functions for jets with pT reconstructed in the range 
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Figure 3: Color online: The reconstructed jet energy reso- 
lution determined from PYTHIA simulations (histograms) 
and real di-jet data for the three jet algorithms used, 
\ri\ <1-R, R=0.7. Red triangles - kr, blue circles - anti-kr 
and SISCone - black squares. 
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Figure 5: Color online: Charged particle, detector level, 
5 fragmentation functions for jets reconstructed with 20< 
PT < 30 GeV/c compared to PYTHIA for 3 different jet 
algorithms. |7;| <1-R, R=0.4. Red triangles - kr, blue 
circles - anti-kr and SISCone - black squares. 



20-30 GeV/c are shown for jet resolution parameters 
R=0.4, Fig.|4]and Fig. [5] and R=0.7, Fig.|6]and Fig.[7| 
Figure |4] and Fig. [6] are the fragmentation functions 
as a function of z {=p^'^'^^°^ /pi^*) and Fig.jsjand Fig.]?] 
as a function of ^(=ln(l/z)). The solid points are the 
data and the histograms are the PYTHIA simulations. 
There is reasonable agreement between the data and 
PYTHIA, and the different jet algorithms (shown as 
different colors and line types/shapes in the figures) 
reconstruct the same fragmentation functions within 
errors. This agreement, especially for the larger res- 
olution parameter, suggests that at RHIC energies 
there are only minor NLO contributions beyond those 
approximated in the PYTHIA LO calculations. 




Figure 4: Color online: Charged particle, detector level, 
z fragmentation functions for jets reconstructed with 20< 
PT < 30 GeV/c compared to PYTHIA for 3 different jet 
algorithms. \ri\ <1-R, R=0.4. Red triangles - kr, blue 
circles - anti-kr and SISCone - black squares. 
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Figure 6: Color online: Charged particle, detector level, 
z fragmentation functions for jets reconstructed with 20< 
PT < 30 GeV/c compared to PYTHIA for 3 different jet 
algorithms. I?;! <1-R, R=0.7. Red triangles - kr, blue 
circles - anti-ky and SISCone - black squares. 



4. The Underlying Event 

The underlying event is defined as everything but 
the hard scattering. Thus, it has contributions from 
soft and semi-hard multiple parton interactions, initial 
and final state radiation and beam-beam remnants, 
Fig. T] Pile-up (multiple events being recorded as one 
by the experiment) is not considered to be part of the 
underlying event. This study is performed at mid- 
rapidity, hence the beam-beam contribution is mini- 
mal. We follow the CDF technique |T3] to perform 
this study. Briefly, the method follows the following 
steps: First the jets in an event are reconstructed, us- 
ing the algorithms mentioned above. Next, the event 
is split into four sections defined by their azimuthal 
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Figure 7: Color online: Charged particle, detector level, 
^ fragmentation functions for jets reconstructed with 20< 
Pt < 30 GeV/c compared to PYTHIA for 3 different jet 
algorithms. \ri\ <1-R, R=0.7. Red triangles - kr, blue 
circles - anti-kr and SISCone - black squares. 
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Figure 8: Cartoon of the areas in A(j>, relative to the lead- 
ing jet, used to define the near jet , away jet , TransMin 
and TransMax regions for the underlying event analysis. 
Figure taken from [5]. 



angle with respect to the leading jet axis (A(/)). The 
range within |A0|<6O° is the lead jet region. An away 
jet area is designated for |A0|>12O'^. This leaves two 
remaining transverse sectors of 60°<A(/)<120° and 
— 12O°<A0<-6O", see Fig. [s] The transverse region 
containing the largest charged particle multiplicity is 
assigned the TransMax region, the other transverse 
sector is termed the TransMin region. The probabil- 
ity of the TransMax region containing contributions 
from initial and final state radiation components of 
the hard scattered parton is enhanced via this selec- 
tion criteria. A "leading" jet study, where at least 
one jet is found in STAR's acceptance, and a "back- 
to-back" jet study, which is a sub-set of the "leading" 
jet collection are performed. This "back-to-back" sub- 
set of events has two (and only two) found jets with 
^^wayjet^^ieadjet^Q^^ and | A0,et | > 150" . This Selec- 
tion suppresses the probability that the scattered par- 
ton undergoes any significant hard initial and final 
state radiation. Thus, by comparing the TransMax 
and TransMin regions in the "leading" and "back-to- 
back" sets we can alter the sensitivity of the data to 
the various components in the underlying event. 

Figure [9] shows the uncorrected charged particle 
densities of the away jet, TransMin and TransMax 
regions of the "back-to-back" data set using the kr 
algorithm, R=0.7. The first observation is that the 
underlying is largely independent of the jet energy 
while the away jet particle density steadily rises with 
reconstructed jet pr, as expected. Figure [TO] shows 
the uncorrected mean px of these regions. Again the 
underlying event results are essentially independent of 
the hard scattering energy. The data for the under- 
lying event at RHIC energies are slightly lower than 
those observed for the underlying event at CDF [T3] 
. The results are consistent within errors for the SIS- 
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Figure 9: Color online: The uncorrected charged particle 
density in the "back-to-back" data set for the away jet, 
the TransMin and the TransMax regions as a function of 
reconstructed lead jet pr, using kr algorithm, R=0.7. The 
dashed histograms indicate the predictions from PYTHIA. 



Cone, kr, and anti-k^ algorithms. In both figures 
the green dashed histogram is the prediction from 
PYTHIA. This PYTHIA tuning gets the trends of the 
data essentially correct, but slightly under predicts the 
particle multiplicity and mean px in the Transverse re- 
gions, while over predicting the activity within the jet. 
In this tuning, as stated above, the multiple parton 
scattering energy scaling factor, e, is set to 0.25. The 
data is still slightly above PYTHIA supporting the ev- 
idence that a large value of the energy dependence pa- 
rameter, e, is needed than that initially determined of 
e=0.16. This also suggests that a number of the LHC 
predictions, often labeled as "ATLAS" tunes, are over 
predicting the activity in the underlying event at 14 
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Figure 10: Color online: The uncorrected charged particle 
mean pr in the "back-to-back" data set for the away jet, 
the TransMin and the TransMax regions as a function of 
reconstructed lead jet pr, using kr algorithm, R=0.7. The 
dashed histograms indicate the predictions from PYTHIA. 
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Figure 11: Color online: The uncorrected charged particle 
density in the TransMin and TransMax regions as a func- 
tion of reconstructed lead jet pr, using SlSCone algorithm, 
R=0.7. 



Figure 1 1 shows the measured charged particle den- 
sity in the underlying event. The densities are the 
same within errors for the "leading" and "back-to- 
back" datasets. This measurement, like those shown 
in Fig. |4]|7j can be used to argue that the hard scat- 
tered partons emit very small amounts of large angle 
initial/final state radiation at RHIC energies. This 
is very different in 1.96 TcV collisions where the 
"leading" / "back-to-back" density ratio is ^0.65 [13]. 
The two solid lines show the expected density if events 
follow a Poisson distribution with an average of 0.36. 
The similarity of this simple simulation to the data 
suggests that at RHIC energies the splitting of the 
measured TransMax and TransMin values is predom- 
inantly due to the sampling. PYTHIA results, not 



shown here, again display satisfactory agreement with 
the data for both data sets. 



5. Summary 

In summary, jet fragmentation functions have been 
measured in p-p collisions at ^/s = 200 GeV and will 
provide a stringent baseline for the measurements un- 
derway in Au-Au collisions. The three jet algorithms 
studied, SISCone, k^ and anti-ky provide consistent 
results. PYTHIA, tuned to 1.8-1.96 TeV data, shows 
reasonable agreement suggesting that the the underly- 
ing physics and its energy dependence is well modeled. 

The underlying event is largely decoupled from 
magnitude of the momentum transfer of hard scatter- 
ing. The data show there is virtually no large angle 
initial and final state radiation from the hard scat- 
tering at RHIC energies. The collision energy depen- 
dence of the multiple parton interactions in the event 
is more significant than initially estimated. This re- 
sults in a smaller prediction for the magnitude of the 
underlying event contribution to the background of 
the reconstructed jets at LHC energies. At RHIC the 
Pt spectra of particles in the underlying event are sig- 
nificantly softer than those from jets. 

In the future these data will be corrected for detec- 
tor inefficiencies and irresolution allowing for compar- 
ison to theoretical models at the particle level. The 
structure of the jet shapes (k^, Jt etc) will be stud- 
ied, as will the particle type composition of the jets 
and the underlying event. Ultimately all these results 
will be compared to similar measurements made in 
heavy-ion collisions at the same beam energy. 
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